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ABSTRACT 


The Unmanned Aerial Vehicle (UAV) plays a critical role in the current 
battlefield in the areas of information superiority, collateral damage, urban area fighting 
and precision strike against high payoff targets. The desire for a shorter “kill-chain” is 
driving the evolution of the Unmanned Combat Aerial Vehicle (UCAV). However, due to 
high cost and limited quantities, UCAVs are currently only available to military planners 
at the operational or strategic level. 

This thesis aims to provide a low-cost solution through integrating commercial 
off-the-shelf (COTS) technologies to produce a prototype of a “Tactical Unmanned 
Combat Aerial Vehicle UCAV” system that can be utilized by the front-line ground units 
in the near future. The Tactical UCAV is designed to enhance the information collection 
and autonomous precision strike capability of the ground units. The Tactical UCAV can 
also be deployed as sensor nodes as part of a larger global information grid in a network- 
centric warfare operation. 

The proposed Tactical UCAV system is comprised of a Hunter Unmanned Aerial 
Vehicle (HUAV), which primarily carries high resolution sensors and communication 
devices and is used as a mother-ship for smaller “Killer UAVs (KUAV).” The KUAV 
carries a mission specific set of instruments; it can be a sensor or a warhead or both 
depending on the desired end results. After the target is acquired by the HUAV, the target 
information will be transferred to the KUAV. The KUAV can then be launched in close 
proximity of the target with the target position update from the HUAV. 

This thesis will focus on the development of a prototype KUAV and the 
integration of the prototype with the existing HUAV “Rascal” developed and operated by 
the Naval Postgraduate School (NPS). The KUAV and the HUAV will form the Tactical 
UCAV system. 
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I. INTRODUCTION 


A. OPERATIONAL REQUIREMENTS 


Combat missions in the past ten years have proved beyond any doubt that the 
dependence on the Unmanned Aerial Vehicle (UAV) has increased significantly. The 
demand for UAVs has evolved from conventional hot-war engagements to low-intensity 
conflicts such as peacekeeping and humanitarian missions. The demand is so great that it 
has outpaced the ability for research and development of UAVs to meet the wide 


spectrum of operational requirements. 


The UAVs in the current battlefield play a critical role in areas such as 
information superiority, collateral damage, urban area fighting and precision strikes 
against high payoff targets. The three main areas that have been identified for future 
UAV developments include the growth in size of strategic UAVs for greater payload and 
endurance, reduction in size of tactical UAVs and the addition of weapons to UAVs to 


offer lethal capabilities in combat missions [1]. 


The desire of shorter “kill-chain”—the time from when the target is spotted by a 
sensor until a shooter locks on—is the main force driving the evolution of the Unmanned 
Combat Aerial Vehicle (UCAV). The current trend towards asymmetrical warfare, where 
90% of the targets are mobile or of the “emerging” variety, emphasizes this even more. 
The effects of a short “kill-chain” are less apparent at the operational and strategic levels, 
while the values of rapid targeting are obvious at the tactical level of war. If the process is 
slow, time critical opportunities may be lost or mobile targets may be long gone by the 
time an air strike arrives. The faster the attack, the less time the enemy has to react, adjust, 


adapt, mount a counter-offensive or escape [2]. 


This thesis aims to provide a low-cost solution through integrating commercial 
off-the-shelf (COTS) technologies to produce a tactical hunter-killer system that satisfies 
the needs for both the reduction in size of tactical UAVs and the weaponization of UAVs 
to offer lethal capability. 


B. PROPOSED CONCEPT OF OPERATIONS 


The proposed tactical hunter-killer system is better known as Tactical Unmanned 
Combat Aerial Vehicle or Tactical UCAV. Currently, due to high cost and limited 
quantities, UCAVs are only made available to military planners at the operational or 
strategic level. Low-cost production through integrating commercial off-the-shelf 
technologies will be the driver for making such superior war-fighting capabilities 


accessible to ground units in the near future. 


The Tactical UCAV is designed to enhance the information collection and 
autonomous precision strike capability of the ground units. The Tactical UCAV can also 
be deployed as sensor nodes as part of a larger global information grid in a network- 
centric warfare operation. The proposed concept of operations for the Tactical UCAV is 


shown in Figure 1 and can be summarized as follows: 


A Hunter Unmanned Aerial Vehicle (HUAV) carries high-resolution sensors and 
communication devices. As the name implies, the HUAV will predominantly play the 
hunter roll in the proposed Tactical UCAV system. She is configurable to carry one or 
many of Killer Unmanned Aerial Vehicles (KUAV), and in this aspect she can be thought 
of as the “Mother Ship” or the air command center for the KUAVs. 


The KUAV can be configured to carry a payload of a sensor, a warhead or both, 
and it will predominantly play the roll of the killer in the proposed system. 


After the target is located by the hunter, the information about the target will be 
“transferred” to the killer. The killer will then be autonomously launched from the mother 
ship. The killer vehicle will autonomously maneuver to the target with the help of 


constant updates from the mother ship as to the position of the target. 
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Figure 1. Proposed concept of operations 


The Tactical UCAV system is designed to be configurable and low cost, so it can 
be deployed as a single system or in numbers (swarm) to fulfill many operation 


requirements. The possibilities include: 
Single System 


Hunt and Kill Mission. As the name implies, upon acquiring a high-payoff target 
by the hunter, a killer will be launched to destroy the target. 


Continuous Surveillance Mission. Instead of destroying the target, the “killer” 





carrying surveillance equipment can be launched to close in or follow the target for 


continuous surveillance. 
Multiple Systems (Swarm) 


Hunt and Kill Mission. This is similar to the first single system mission. 





However, it would be used when multiple targets, such as a column of tanks, are to be 
engaged simultaneously. Multiple Tactical UCAVs can be assembled in the area of 
operation. If each launches two killer vehicles, a swarm of killers is now available to 


simultaneously engage the targets. 


Continuous Surveillance Mission. A single system has a limited area of coverage 





in any surveillance mission. A swarm of sensors forming a sensor grid can be deployed to 
cover a large area or track a moving target with greater tracking accuracy. This operation 
can also be part of the global information or sensor grids in a network-centric warfare 


operation. 
C. SCOPE AND OBJECTIVES 


It can be seen from Figure 1 that this thesis will focus on the development of a 
prototype KUAV. Due to the limited duration available for this thesis, the integration of 
the prototype KUAV with the existing HUAV “Rascal” [3] developed and operated by 
the Naval Postgraduate School (NPS) will be done in the near future. The combination of 
the prototype KUAV with the existing HUAV will form the Tactical UCAV system. 


The scope of this thesis includes the following: 

Development of Prototype 

Design and build the test prototype 

Aerodynamic performance analysis with linear panel software 
Tradeoff studies for various performance parameters 
Guidance and control analysis with numerical methods 
Definition of engagement envelope 

Hardware-in-the-Loop Testing 

Architecture and setup 

Procedures for working with Procerus Virtual Cockpit 2.4.2 


Refining engagement envelope 
D. THESIS ORGANIZATION 


Chapter II provides background on the UCAV and Tactical UAV systems 


currently in service or development. Chapter II presents the critical considerations and 
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decisions made for the design and building of the prototype. In Chapter IV, the 
methodology and tools used for generating aerodynamics performance data from linear 
panel codes will be discussed. Subsequently, the results of the tradeoff studies for key 
performance parameters of the prototype will be presented. In Chapter V, the 
methodology and simulations used for modeling the practical operating environment with 
numerical techniques will be discussed. The results of the simulations will also be 
presented in this chapter. In Chapter VI, the setup, procedures, results and findings 
obtained from Hardware-in-the-Loop (HIL) experiments will be discussed. Finally in 


Chapter VII, the conclusion and recommendations for future works will be discussed. 
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I. LITERATURE REVIEW 


A. EXISTING UCAV SYSTEMS 


Major UAV systems have been well developed over the years. In recent years, 
much work and effort have been invested in the development of hybrid systems to reduce 
the time required in the sensor-to-shooter targeting cycle. These hybrid systems primarily 
integrate the sensor and shooter in a single platform more commonly known as 
Unmanned Combat Air Vehicles (UCAV). As part of the literature review, two in-service 


UCAVs are presented. 
1. Predator B — General Atomics Aeronautical Systems, USA 


The Predator B aircraft was developed by General Atomics Aeronautical Systems, 
USA, in 2000 and its flight commencing in February 2001. Powered by a turboprop 
engine, the Predator B series was designed as a long endurance, high-altitude unmanned 
aircraft for use as a multi-mission system. From reconnaissance, surveillance, targeting 
and weapons delivery to scientific research and other civilian applications, Predator B has 
the capacity to conduct multiple missions simultaneously due to its large internal and 
external payload capacity. The detailed technical specifications are presented in 
Appendix A. [4] 





Figure 2. Predator B armed with Hellfire missiles 
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2. Sperwer B — Sagem, France 


Sperwer B was designed and developed by Sagem, France, to support Intelligence, 
Surveillance, Target Acquisition and Reconnaissance (ISTAR) at the battlegroup level. 
Sperwer is capable of carrying multiple payloads ranging from Forward Looking Infra- 
Red (FLIR), Electro Optic / Infra-Red (EO/IR) and Synthetic Aperture Radar (SAR). It is 
equipped with two underwing hardpoints that can carry external loads of up to 30kg each. 
Sagem has already demonstrated the integration of Sperwer B with the Spike LR missile. 
They are cooperating with GIAT, to test a new smart munition delivery system, based on 
the Bonus submunition. However, the armed configuration requires that up to 20 kg of 
fuel must be removed. This will limit the endurance of the armed Sperwer. The detailed 


technical specifications are presented in Appendix B. [5] 





Figure 3. Sperwer B armed with Spike LR missiles 


B. EXISTING TACTICAL UAV SYSTEMS 


The development of the Tactical UAV has also progressed significantly in recent 
years. They are evolving into multi-role, multi-mission platforms. As UAV technology 
matures, UAV will become increasingly cost effective, they will grow smaller and be 
able to accomplish a greater number of missions and fill a greater number of roles. 
Besides their current applications in Reconnaissance, Surveillance and Target Acquisition 


(RSTA), the Tactical UAV mission set could be expanded to include target designation, 


strike, mine countermeasures, electronic warfare and information warfare. Also as part of 


the literature review, two in-service tactical UCAVs will be presented. 
if Shadow 200 —- AAI Corporation, USA 


The Shadow 200 short-range Tactical UAV was developed by AAI Corp. in the 
1990s. It has a non-retractable tricycle landing gear for conventional wheeled take-off 
and landing. The RQ-7A can also be launched from a catapult and has a tailhook to catch 


arresting cables for a shorter landing strip. 


Shadow 200 is capable of spotting targets up to 125 kilometers away from the 
brigade tactical operations center. It has both day and night capability of identifying 
tactical vehicles up to 8,000 feet above ground at a slant range of more than 3.5 
kilometers. The Shadow ground control station provides near real-time targeting data for 
precision weapons using “leap ahead” technology and transmits near real time imagery 
and telemetry data to the end users. The detailed technical specifications are presented in 


Appendix C. [6] 





Figure 4. | Shadow 200 being launched from a catapult 


2; Watchkeeper Tactical UAV — Thales, UK 


The Watchkeeper Tactical UAV system will provide the UK armed forces with 
Intelligence, Surveillance, Target Acquisition and Reconnaissance (ISTAR) capability. 
Watchkeeper is a tactical system that will be operated on the battlefield by the British 
Army Royal Artillery. The air vehicle will be capable of carrying a range of sensors 
including day and night cameras and surveillance radars. Two WK450 air vehicles will 
be able to operate in tandem, with the second acting as a communications relay. The 
ground control station will be network enabled to ensure comprehensive communications 
links to, for example, airborne stand-off radar, attack aircraft and battlegroup 
headquarters. The Watchkeeper system will enter service in the British Armed Forces 
Royal Artillery in 2010. The detailed technical specifications are presented in Appendix 
D. [7] 
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Figure 5. | Watchkeeper Tactical UAV 
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HiIl. BUILDING THE PROTOTYPE 


A. PROTOTYPE DESIGN 


After much research and many reviews have been done on open-source data, it is 
apparent that there is no commercially available system that can meet the requirements of 
the proposed KUAV. Therefore, part of this thesis involved the designing and building of 


the prototype that is capable of integrating with “Rascal” for concept demonstration. 


The proposed concept of operations placed many constraints on the design of the 


KUAV. The three principal considerations for the KUAV are as follows: 


Size and weight—It must be small and light so that it can easily integrate with 


“Rascal” 


Lift and drag—It must have high lift and low drag in order to increase the payload 


and reduce the thrust requirements 


Noise signature—The noise signature must be kept to a minimum so that it can 


have a stealthy approach 


With the principal considerations in mind, a non-powered glider with swept 
linearly tapered wing has been selected as the design for research and experimentation in 


this thesis. 


Due to the tight time schedule, a quick study was conducted on in-service glider 
systems currently available. The decision was made to design the prototype by scaling 
down (maintaining an aspect ratio of AR = 2.78) the “Swift 2 Wing” hobby flying wing 
manufactured by Northeast Sailplane Products [8]. Apart from being able to rapidly 
produce the prototype, the selection of the autopilot was the other important factor. 
KESTREL™ Autopilot 2.2, manufactured by Procerus Technologies, [9] will be used to 


control the autonomous flight of the KUAV. This autopilot is being used on a similar 
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class of test platforms that have tapered wings that are swept linearly. Therefore, this 
could potentially save much effort and time on tuning the Proportional, Integral and 


Derivative (PID) gains of the autopilot. 
1; Dimensions 


Apart from the reasons cited above, the size of the prototype also ensures that it 
can be easily carried and released without hindering the normal operations of “Rascal.” 
EH 3.0/12 (3% camber and 12% thick) airfoils, designed by [10], were selected. Table | 
shows the detailed dimensions. Figure 6 shows the plan view of the prototype and the 


cross section of the airfoil. 


Root Chord Tip Chord 


ip 
| (mm) | (mm) | (%) | (mm)_| (deg) |_ (mm) _| 


Table 1. | Dimensions of prototype 





Figure 6. _ Plan view of prototype and cross section of EH 3.0/12 
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2. Material 


The wing was made from commercial off-the-shelf hobby grade Spyder foam. 
The foam is durable, lightweight and more importantly, readily available at a low cost. It 
is can be easily cut and shaped using a heated wire. The two main disadvantages of 
cutting this material (by amateurs) are the surface roughness and the difficulty of shaping 
a consistent and accurate airfoil shape. They both will decrease lift and increase parasite 
drag on the airfoil, thus reducing the performance of the airfoil. With these disadvantages 
in mind, this is still the best solution for rapid prototyping. The elevons of the prototype 


were made of lightweight balsa wood and the winglets were made of 1/32” aircraft ply. 
3. Components and Loading 


The prototype was first tested for its aerodynamic performance using a Remote 
Controlled (RC) flight mode. Only after successful RC flights was the prototype tested 
for autonomous flight using the Auto Pilot (AP). Different commercial off-the-shelf, 
hobby-grade components were required for each of these individual flight modes. A 
detailed breakdown of the components and the respective weight for both modes are 


presented in Appendix E. 
4. Trimming and Balancing 


The theoretical aerodynamic center of the prototype was estimated to be at 25% 


of the Mean Aerodynamic Chord (MAC) computed using equation (1). 


1 2 
ieee) 
3 (1412) 


where rc=Root Chord and t=Taper Ratio 


(1) 


MAC was computed to be 181 mm and the theoretical aerodynamic center is 
located at 128 mm on the root chord from the nose. After balancing all the components 


and trimming from hand-launched RC test flights, the optimum center of gravity CG was 
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found to be located at approximately 118 mm from the nose on the root chord resulting in 


a static margin of 10 mm. The details are depicted in Figure 7. 








25% of MAC 






Center of Gravity --~~ 


Aerodynamic Center ~~~ 


Figure 7. Location of aerodynamic center and CG 


Since this static margin is already sufficient for trim flight during hand-launched 
test flights (low altitude and short duration). The author is confident that the prototype 
KUAV will exhibit more stable performances during “Rascal” launched test flights (high 


altitude and long duration). 
5. Completed Prototype 


The main components (from left to right: autopilot, battery-pack, GPS and 


modem) used for the autonomous flight are shown in Figure 8. 





Figure 8. | Components for autonomous flight 
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The prototype KUAV is shown Figure 9. The components described above except 
the GPS, will be packed in an opening cut on the airfoil under the plastic cover shown in 


Figure 9. A separate slot will be cut on the airfoil to house the GPS unit. 





Figure 9. Side and top views of the prototype KUAV 
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IV. AERODYNAMICS STUDIES 


A. BASIC AERODYNAMICS FOR GLIDERS 


The aerodynamic forces acting on a non-powered glider closely resemble those on 
any aircraft. The key distinction is that the non-powered glider does not produce thrust 
force to counteract the drag force. As a result, a glider cannot maintain a level flight at 
constant speed and altitude. A glider in steady gliding flight is always descending relative 
to the air around it. The glider is essentially trading altitude to maintain its velocity. 
Figure 10 shows the balance of forces acting on a glider descending in a steady glide at 


constant velocity. 


Lift (L) 


Horizontal 


Glide Angle (y) ( “a 


- 
- 


Drag (D) _--7 


-_<-- 
ee 





---c7 
ia 
-- 







Angle of Attack (a) 


Velocity (V) — Flight Path 


Figure 10. Forces acting on a glider in flight 


The glider follows a path that slopes at an angle below horizontal called the Glide 
Angle, denoted by the Greek letter (y). The direction of the glider’s motion along the 
glide path is shown by the velocity arrow on the diagram. (Note that the glider’s 
longitudinal axis is not pointing along the direction of motion, but slightly above the 


glide path at angle a, called the Angle of Attack, AOA). 
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The lift force (LZ), which according to linear theory acts at a right angle to the 
velocity through the air, is tilted forward from vertical at an angle equal to y. Drag (D) 


acts directly opposite the direction of motion through the air. 


At trim condition (constant glide angle and velocity) the following relationships 


can be derived. 
L=W, — L=mgcosy (2) 
D=W, — D=mgsiny (3) 


It can be observed that when all forces on the glider are balanced there is no net 
force to cause acceleration, so the glider will move along the glide slope at a constant 
velocity. The velocity (V) and glide angle (vy) can be derived from equations (2) and (3) 


respectively. 


Le = PAV*C, and D= = PAV'C, 








: pAV’C, =mg cosy (4) 

gees . 

aroY C, =mgsiny (5) 
tany=—2 +> y=tan"' i > y=tan" = (6) 

L L L 
2 
V= —_— (7) 
pA(C, cosy +C, siny) 


B. SIMULATION USING ANALYTICAL TECHNIQUES 


Commercial off-the-shelf software, LinAir (Version 4.3) [11] was used in this 
thesis to analyze the aerodynamic characteristics of the glider. LinAir is a linear panel 
code program capable of computing aerodynamic characteristics of multi-element, 


nonplanar lifting surfaces for a given wing geometry and angle of incidence. The results 
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generated are based on solving the Prandtl Glauert linear partial differential 


equations for inviscid, irrotational and subsonic flow. 
1. Verification of Panel Codes 


There is a need to first “calibrate” several input parameters to LinAir, so that the 
results generated can be benchmarked with existing wind tunnel results for a similar class 
of wing profile. This is especially true for cambered airfoil, as the procedures for 
modeling airfoil with camber stated in [11] are strictly adhered to. It was found that 
EH3.0/12 has zero lift at AOA of approximately -1°. Therefore, all results generated will 
be adjusted by adding -1° to the actual incidence angle. 


In the proposed operating envelope of the glider (low altitude and low speed) the 
Reynolds Number (Re) is expected to fall in the region of 50,000 to 100,000. Wind tunnel 
experimental results obtained by [12] provide a good insight into the aerodynamic 
performance of low aspect ratio wings at low Re. This served as an excellent independent 


source to cross-check the data generated by LinAir. 


Many aerodynamic parameters, in particular those related to drag, cannot be 
computed analytically or by linear methods and are even difficult to get right using 
Navier-Stokes solvers, and are usually obtained from wind tunnel experiments. Extensive 
tests for wings of AR = | and 2 at Re = 70,000 and 100,000 have been conducted and the 
results are presented in [12]. In the absence of wind tunnel facilities and time, the 
following parameters, Drag Coefficient at Zero Lift (C,, ) and Induced Drag Coefficient 
(K ) for the prototype glider (AR = 2.78) will be inferred based on the results published 
by [12]. 

It is assumed that the total drag produced by a wing is give by equation (8). The 
values of C, and K can be obtained by plotting C,, vs. C,” for each model and applying 


a least-square linear regression to the data. The wind tunnel data obtained by [12] are 


summarized in Table 2. 


Cy =Cy, + KC, (8) 
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1 0.015 
2 0.020 


Table 2. Wind tunnel results of C,, and K for wings of AR | and 2 





These data were subsequently run on LinAir using wing profiles matching the 
respective AR. The results generated for AR = 2 (AOA from 0° to 10°) closely matched 
those obtained in the wind tunnel. However, for AR = 1, LinAir produced a better match 


using K =0.4. The detailed results are presented in Appendix F. 

Since the prototype has an AR of 2.78, by extrapolating from the above simple 
verification, C,, = 0.022 and K =0.4 will be used to model the prototype glider. LinAir 
has proven to be suitable and will be used as the analytical tool for estimating the 
aerodynamic performance of the prototype. 


2. Geometry and Input File 


The geometry of the glider described in Chapter III was modeled using a total of 
AO panels covering the entire wing surface and the two elevons as shown in Figure 11. 


The input file is shown in Appendix G. 














Figure 11. Geometry of glider generated by LinAir 
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3. Aerodynamic Data 


The aerodynamic performance of the glider was evaluated for AOA from 0° to 20° 
at an interval of 0.5°. (This range of AOA is only feasible for ideal linear theory 
simulations using LinAir. The operating AOA range is expected to be lower in practical 
model and physical implementation.) Negative AOA were not examined as it is not 
possible to maintain steady flight for negative AOA for a non-powered glider. More 
importantly, it does not contribute to any performance improvement. Numerical results 
generated by LinAir for both the non-cambered airfoil and with the adjustments made for 


the cambered airfoil are shown in Appendix H. 
C. ANALYTICAL TRADEOFF STUDIES 


Three main performance parameters of the glider were examined. 
Glide Velocity (V) 
Horizontal Range (R) 


Airborne Endurance (E£) 


These performance parameters are also functions of AOA and y. Using the data 
generated, the following five tradeoff studies were analyzed at a release altitude of 300 m 


which is a typical operating ceiling for tactical systems. 
Tradeoff Study 1 - L/D vs. AOA 
Tradeoff Study 2 — y vs. L/D 
Tradeoff Study 3 — R vs. L/D 
Tradeoff Study 4— V vs. AOA 


Tradeoff Study 5 — E vs. V 
1. Tradeoff Study 1 — L/D vs. AOA 


It is important to recall at this juncture that the values of L used in this analysis 


are estimated by LinAir. LinAir is a linear panel code program, it assumes that L always 
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varies linearly for all AOA and does not include the effects of flow separation and stall 
which are dominant at higher AOA. Conversely, the effects of flow separation and stall 
are approximated in D, as D is a quadratic curve fit to empirical data. (The plot showed in 
Figure 12 will change if L is computed accurately with considerations to the effects of 


flow separation and stall.) 


This study is confined to the region of AOA around 3° to 6° where the effects of 
flow separation and stall are negligible. Therefore, the results presented in Figure 12 in 


the region of AOA around 3° to 6° are valid to this study. 


Tradeoff Between L/D and AOA 
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Figure 12. Tradeoff between L/D and AOA 


Analysis of Figure 12 shows that the maximum L/D of approximately 4.4 
occurred at an AOA of 3°. The data presented in Table 3 shows the comparison between 


the analytical results to those obtained by [12] through wind tunnel experiments. 


22 


Experimental Wind Tunel Data 
Be New ee ee ere arene 
| 1 | 1.00 | 48 | 6.0° | 70,000 


70,000 
Analytical LinAir Result 
70,000 





Table 3. | Comparison between experimental and analytical results 


The above comparison shows that the results obtained from LinAir are 


comparable to the wind tunnel experiment and the selected values of C, =0.022 and 


K =0.4 are within acceptable engineering estimates. 
2. Tradeoff Study 2 - y vs. L/D 


It was shown in Equation (6) that y is a function of L/D. The greater the L/D, the 


gentler the glide slope will be and vice versa. This result is shown in Figure 13. 
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Figure 13. Tradeoff between y, L/D and AOA 
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The graph of Tradeoff Study 1 was projected on the x-y plane and the z-axis 
shows the values of y with varying AOA and L/D. 


Therefore, following the results of Tradeoff Study 1 and the relationship 
presented in equation (6), y would first decrease with increasing AOA. After stall 
condition, y would increase with increasing AOA. A minimum glide angle of 12.7° would 


occur at the AOA of 3° where L/D was at the maximum. 


3. Tradeoff Study 3 — R vs. L/D 


Horizontal 





Velocity (V) 


R 


Figure 14. Glide path geometry 


The glide path shown in Figure 14 can be related through equation (9). The 
gentler the glide slope (smaller y), the farther the glider will travel horizontally before 
reaching the ground. Therefore, maximum R will occur at minimum y. 


_ A 
tan v 





(9) 


ae > R 
R 
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Figure 15. Tradeoff between R, L/D and AOA 


The result of maximum R occurring at minimum y is shown in Figure 15. The 
graph of Tradeoff Study 1 was again projected on the x-y plane and the z-axis shows the 
values of R with varying AOA and L/D. 


Therefore, following the results of Tradeoff Study 1 and 2 and the relationship 
presented in equation (9), R would first increase with increasing AOA. After stall 
condition, y would decrease with increasing AOA. A maximum R of approximately 1.335 


km would occur at the AOA of 3° where y is at the minimum and L/D is at the maximum. 


4. Tradeoff Study 4 -— V vs. AOA 


It has been explained that a non-powered glider does not produce thrust to counter 
drag. In order to maintain a trimmed glide path, it is always trading off altitude to 
maintain its velocity. Equation (7) shows the analytical expression of the glide velocity at 


trimmed flight with a fixed AOA. 
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Tradeoff Between V and AOA 
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Figure 16. Tradeoff between V and AOA 


Analysis of Figure 16 shows that V increases exponentially at low AOA and 
decreases gradually at high AOA. This behavior of V is critical to the next tradeoff study 


where speed is varied to achieve the desired airborne endurance. 
5, Tradeoff Study 5 - E vs. V 


Trimming the glider to fly at its maximum L/D will allow it to fly farthest 
horizontally, but it will not result in the longest airborne endurance (£) (duration of time 


in the air). 


E depends on the vertical sink rate, or rate of decrease in altitude (H), as shown in 
equation (10). This depends on both the glide angle and the glide speed. Maximum E is 
achieved by trimming the glider to fly at an AOA higher than that where maximum L/D 


occurs. (L/D will decrease slightly from the maximum value.) This will result in a steeper 
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glide slope, but it will cause a significant reduction of glide velocity. The reduction in 
speed is sufficient to compensate for the effects of the steeper glide slope that are needed 


to achieve maximum EF. 


H . 
a where V,=Vsiny (10) 
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Figure 17. Tradeoff between E, V and AOA 


Maximum FE is achieved by trimming the glider to fly at a higher AOA, and this 
result is shown in Figure 17. The graph of Tradeoff Study 4 was projected on the x-y 
plane, and the z-axis shows the values of E with varying AOA and V. 


Therefore, following the results of Tradeoff Study 4 and the relationship 
presented in equation (10), E would first increase with increasing AOA until slightly after 


maximum L/D condition. Following maximum L/D condition, E would decrease with 
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increasing AOA because a decrease in speed is insufficient to compensate for the effects 
of the increasing glide angle. A maximum E of approximately 124 s would occur at the 


AOA of 6°. 
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Figure 18. Comparison between R and E vs. AOA 


Finally, it can be concluded from the analysis of Figure 18 that maximum R 
occurs at an AOA of 3° where L/D is the maximum. Maximum E occurs at a slightly 
higher AOA of 6° where further reduction in speed is sufficient to compensate for the 
effects of the steeper glide slope. Maximum E can only be achieved at the expense of 


horizontal distance (approximately 16.5% reduction in R). 


The results obtained from the above theoretical tradeoff studies under ideal 
operating conditions will define the most optimistic operational envelope of the prototype 
KUAV. The tradeoff between critical performance parameters provides a basis for the 
end-user to plan for different engagement scenarios like long range (where maximum 


range is desired), loiter (where high endurance is desired) and dash (where high speed is 
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desired). In Chapter V, the operational envelope of the prototype KUAV will be 
examined under practical operating conditions. These conditions will be modeled as 


closely as possible using numerical methods (MATLAB Simulink). 


29 


THIS PAGE INTENTIONALLY LEFT BLANK 


30 


V. GUIDANCE AND CONTROL 


A. SIMULATION USING NUMERICAL TECHNIQUES 


The analysis of the aerodynamic performance data generated by LinAir was based 
purely on analytical closed-form solutions derived from equations of motions. These 
ideal conditions can never be replicated in practical implementations. The actual 
performance of the prototype will be dependent on many other factors such as sensor 
measurement errors, variations in atmospheric conditions, the performance of the 


autopilot and control algorithm and the efficiency of the mechanical devices onboard. 


It is almost impossible to simulate the exact operating environment with current 
technology, in particular for conditions involving low Re where variations in 
aerodynamic behaviors and performances are highly non-linear. However, there is still an 
absolute need to account for as many imperfections as possible in order to better estimate 
and predict the operational envelope of the prototype KUAV. MATLAB Simulink and 
the aerodynamic performance data generated by LinAir will be utilized to model a more 


realistic operating environment for the prototype. 


This process is summarized in Figure 19. 





Analytical Performance in LinAir 


6DOF + AP Model in Simulink 


Modeling Prototype in Simulink 
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Figure 19. Flow chart depicting simulation process 
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1. Analytical Results from LinAir 


Chapter III illustrated how the results generated by LinAir were used to 
investigate the aerodynamic performance of the prototype based on analytical equations 
of motions in ideal operating conditions. These dimensionless coefficients from LinAir 
will be used to define the aerodynamic performance of the prototype in the model 


described below. 


2. 6 Degree of Freedom (DOF) + Auto Pilot (AP) Simulink Model 


An existing Simulink Model for Silver Fox UAU previously developed by [13] 
has been used to model the performance of the prototype. A straightforward modification 
was made to the model to more closely approximate the configuration of the prototype 
KUAV. The engine was removed from the model because the non-powered glider does 
not produce thrust. The modified model is shown in Figure 20. It is comprised of the 


following three main blocks: the AP Inner loop, Forces and Moments and the 6DOF 


Rho WYZ 
ISA Atmosphere 
XYZ TP) Sion 


equations of motion. 
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Figure 20. Simulink Block Diagram for DOF + AP Model 


The 6DOF EOM block implements the equations for a constant-mass rigid-body 
6DOF aircraft with respect to the body frame {B}. The detailed derivations will not be 
presented in this thesis as they can be found in [14]. In summary, the motions of a 6DOF 


aircraft can be completely described by the following set of equations: 
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Linear Momentum Equation 

wa =m( "0+ "Ox *) 

Angular Momentum Equation 

°G=J*a+"ax(J*0) 

Euler Equations relating body rates and attitudes angles 
p= p+q(sing+rcos¢)tané 

6 =qcos¢—-rsing 

y =secO(qsing+rcos0) 


The Forces and Moments block implements the forces and moments acting on the 


prototype due to aerodynamics and gravitational effects. They can be expressed as: 


°F = ?F oro ay a 
°G ‘ Gs: 
For simplicity, the aerodynamic forces and moments acting on the aircraft are 


often defined in terms of dimensionless aerodynamic coefficients that are functions of the 


deflections of the control surfaces. 


Cp 
wR| Cy 








Where gq is the dynamic pressure and S is the wing reference area. 


ba R is the transformation matrix that rotates a vector from the wind coordinate 


frame to the aircraft body frame. 
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The numerical values for the dimensionless aerodynamics coefficients were 
determined from the data generated by LinAir. The detailed explanation will be presented 


in the next section. 
The gravitational forces acting on the rigid body will not generate moments as 


they are assumed to be acting on the center of gravity. 


0 
Bp B 
F grav = ik 0 
mg 
WR is the rotation matrix that rotates a vector from the earth frame to the aircraft 
body frame. 


In summary, the forces and moments acting on a 6DOF aircraft can be completely 


described by the following equations. 





Ce 0 
°F =qS ~R|\ C, |+7R| 0 
Cc. mg 
Cb 
°G = S ,R\ C,,c 
Cb 


Finally, the AP Inner Loop block implements the inner loop of the autopilot 
developed by [13]. 


3. Modeling the Prototype KUAV in the 6DOF + AP Model 


There is a need to initialize the (DOF + AP model at the operating conditions and 
performance of the prototype KUAV. This initialization process is critical so that the 
model will generate accurate results based on realistic operating conditions and 
aerodynamic performance of the prototype KUAV. As mentioned in the previous section, 
the dimensionless coefficients obtained from LinAir will be used to model the 


aerodynamic performance of the prototype. 
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The initialization script file of the model is shown in Appendix I. The coefficients 
of performance under the “Aerodynamic Derivatives” are initialized using the data 
generated by LinAir (units of per radian). The following approximating relationships are 


used to compute the coefficients: 
Cy =Cyp t AC, + AC, 
Cy =Cy,B+C,5,0r 
Cid 


B= 
i 


C, =C,,+C,,@+C,,a4+C,,0e+ 





C, =C,,8+C5,5a+ Cod Taal ppt Ce) 





; Cong IC 
C,, =C yo t+ Cp ht C,H +C,,5,0€ += 


moe B— 
Pr 


b 
C,=C,,B+ CcBet Cu 51+ Ta Cw p+C,,r) 
Where 6e, 6a, 6r represent the deflection of the elevator, aileron and rudder respectively. 
(As the prototype KUAV is designed to be a flying wing, there is no rudder and the 
elevons function as both elevators and ailerons. 6r=0 and de=0a.) 
p.q.r are the angular rates, b is the wing span, cis the mean aerodynamic chord 


and @ and £ represent the angle of attack and the sideslip angle respectively. 


4. Simulation Procedures 


In all realistic non-linear models, huge amount of transient response is often 
expected which is detrimental to the accuracy of the results generated. It is important to 
note at this point, that only the results from the steady-state conditions, better known as 
trim conditions, will be used for analysis. Therefore, there is a need to minimize the 


duration of this transient response. This is achieved by running the simulation with the 
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initial conditions as close to the trim conditions as possible. This will greatly improve the 


accuracy of the results obtained from the simulations. 


Specifically for this study, it has already been shown that the performance of the 
non-powered glider depends on the trim conditions. This is primarily determined by the 
glider’s AOA which in turn is controlled by the deflection of the elevators. Therefore, it is 
necessary to ensure that the initial conditions, with respect to the body frame, correspond 


to trim conditions. The following procedures were adopted for this study: 


The Linearization Task function in Simulink Control Design Blockset was used to 
obtain the steady-state velocity for every 1° change in elevator deflection. It was found 
in this step that steady-state conditions could only be achieved for AOA between 1° and 
8°. This operating range of AOA is significantly smaller than those used in Chapter IV. 
This is expected as LinAir is programmed based on linear theory (linearized solution does 


not diverge) and it will generate results even for an AOA of 90°. 


The body velocity in the initialization script file was set to the corresponding trim 
conditions of the respective AOA. (This is to minimize the transient response and improve 


the accuracy of the simulation results.) 
The model was run for AOA from 1° to 8°. 


The output from the workspace for Range (x-position), Velocity and Endurance 


(time) was saved when the altitude (z-position) reached zero. 


The results of the simulations for a release altitude of 300m are attached in 


Appendix J. 
5; Results and Comparisons 


The results obtained from the 6DOF + AP model were compared to the results 


found in the ideal operating environment (presented in Chapter III). 
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Comparison Between Analytical and (DOF + AP Models for L/D vs. AOA 
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Figure 21. Comparison between analytical and 6DOF + AP models for L/D vs. AOA 


Analysis of Figure 21 shows that the L/D of the 6DOF + AP model is about 16% 
lower than the L/D obtained in the analytical model. This result was expected, as 


imperfections in practical environments will reduce the performance of any system. 
Maximum L/D = 4.4 at AOA = 3.0° (Analytical model) 


Maximum L/D = 3.7 at AOA = 2.8° (6DOF + AP model) 
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Comparison Between Analytical and 6DOF + AP Models for V vs. AOA 
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Figure 22. Comparison between analytical and 6DOF + AP models for V vs. AOA 


Analysis of Figure 22 shows that the change in V with increasing AOA exhibits 
the same behavior in the 6DOF + AP model as the analytical model. In the region 
(AOA<2.5°) where V decreases exponentially, system imperfections caused the glider to 
travel at a higher speed. In the region (AOA>2.5°) where V decreases gradually, system 


imperfections caused the glider to travel at a lower speed. 
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Comparison Between Analytical and (DOF + AP Models for R vs. AOA 
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Figure 23. Comparison between analytical and 6DOF + AP models for R vs. AOA 


Analysis of Figure 23 shows that the R generated from the 6DOF + AP model is 


about 17% lower than the R obtained from the analytical model. 
Maximum R = 1.335 km at AOA = 3.0° (Analytical model) 


Maximum R = 1.110 m at AOA = 2.8° (6DOF + AP model) 
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Comparison Between Analytical and 6DOF + AP Models for E vs. AOA 
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Figure 24. Comparison between analytical and 6DOF + AP models for E vs. AOA 


Analysis of Figure 24 shows that the practical E is about 10% lower than the E 


obtained in an ideal environment. 


Maximum EF = 124 s at AOA = 6.0° (Analytical model) 
Maximum E'= 111.1 s at AOA = 4.2° (6DOF + AP model) 
6. Refinement of Operational Envelope 


It was mentioned in Chapter IV that the results obtained from LinAir predict the 
most optimistic (upper bound) operational envelope of the prototype KUAV. The results 
obtained from the 6DOF + AP model in this simulation predict the most conservative 


(lower bound) operational envelope. 
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The maximum horizontal range where the prototype can travel (with a release 
altitude of 300 m) is 1.11 km at an AOA of 2.8°. Alternatively, it can sacrifice range to 


achieve a maximum airborne endurance of 111 s atan AOA of 4.2°. 


The actual performance of the prototype KUAV should fall within the upper and 
lower bound predicted by the analytical and the DOF + AP model respectively. Prior to 
the physical test flight, it is of great importance to conduct the hardware-in-the-loop 
(HIL) experiments to predict the actual performance of the prototype KUAV. The HIL 


procedures and results are presented in Chapter VI. 
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VI. HARDWARE-IN-THE-LOOP EXPERIMENTS 


A. ARCHITECTURE AND SETUP 


Due to the tight timeline of this research, there is no means for the author to 
develop a new guidance law to provide reference commands to the autopilot inner loop. 
This is similar to the AP Inner Loop block described in Chapter V. Therefore, COTS 
Kestrel Autopilot is selected as the hardware for autonomous flight. Its default GPS 
waypoint navigation algorithm will be utilized to provide guidance reference commands 


to the autopilot inner loop. 


The next step to perform is hardware-in-the-loop (HIL) experiments using the 
physical autopilot hardware. These will verify the results obtained in analytical and 
numerical analysis done prior to actual test flight. Kestrel Autopilot and Virtual Cockpit 
have a built-in ability to stmulate a 6DOF UAV through the use of a 3rd party, open- 
source simulator called Aviones. Aviones displays the simulated flight in 3D space, 


which allows the user to quickly and easily visualize the flight plans. 


This is the first attempt by the team to set up HIL for Kestrel Autopilot. Therefore 
it is important that the system architecture, set up and experimentation procedures are 


properly documented for knowledge retention. 
1. Architecture 


There are two options available for setting up the HIL. A brief summary of each 


option is presented. The detailed descriptions are found in the Kestrel Users Guide [15]. 


Architecture 1, which uses a TCP/IP connection as shown in Figure 25. 
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Figure 25. HIL architecture using a TCP/IP connections 
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In the TCP/IP architecture, Virtual Cockpit communicates with Aviones over a 
TCP/IP connection. Aviones passes data (body rates and body velocities) from Virtual 
Cockpit to the autopilot over a hardwired serial connection. The practical advantage of 
this method is that only one serial port is needed (which is the configuration for most 
laptop in the current market). The disadvantage is that the modem needs to be unplugged 
from the autopilot and the results obtained will not account for any latency due to data 


transmission by the modem. 


Architecture 2 which uses a standard serial modem interface as shown in Figure 
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Figure 26. HIL architecture using a standard serial modem interface 


In the serial modem architecture, the autopilot communicates with Virtual 
Cockpit over the standard serial modem interface. Aviones communicates with the 
autopilot over a serial connection. The advantage of this method is that the modem will 


be left connected to the autopilot as per actual operating configuration. The results 
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obtained will give a more accurate representation of the predicted performance. The main 
disadvantage of this set up is that the user will be exposed to unnecessary radiation 


emitted from the commbox. 


Weighing between the advantages and disadvantages of the two options and 


hardware availability, Architecture 1 was selected for this research. 
2. Setup 


In this thesis, the TCP/IP architecture was setup for the HIL experiments. The 


details procedures are presented in Appendix K. 
Additional lessons learned by the author are also documented for future reference. 


When future studies are done, the most current version of Virtual Cockpit 
Aviones must be loaded. (At the time of this experiment, the latest software versions 


were Virtual Cockpit, version 2.4.2 and Aviones, version 2007-11-1) 


When using Virtual Cockpit, version 2.4.2, an update for Microsoft DirectX must 


be installed. The update can be downloaded from 
http://www.microsoft.com/windows/directx/default.mspx 


The most current version of the firmware for the autopilot must be programmed 
into the autopilot hardware. (At the time of this experiment, the firmware version was 


MA8 2 3.6) 
B. | TUNING AUTOPILOT PID GAINS 


Chapter II mentioned that the Kestrel autopilot and Virtual Cockpit were 
developed by Procerus Technologies for the Zagi UAV test platforms. As the prototype 
KUAV is much smaller with different aerodynamic performance than the Zagi, there is a 
need to retune the PID gains of the autopilot to match the dynamic response of the 


prototype KUAV. 
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1. Aerodynamic Performance Parameters 


The dimensionless coefficients from LinAir in Chapter V will again be used to 
define the aerodynamic performance of the prototype in the flight simulator, Aviones. 


The detailed procedures are presented in Appendix L. 
pa Autopilot PID Gains 


After setting up the model in Aviones, the PID gains of autopilot were tuned 
using the HIL simulation. The Ziegler Nichols rules were applied during tuning, with 
addition advice from Chapter II of the Kestrel Installation and Configuration Guide [16]. 
In addition, the author had also found the guidelines provided by Cloud Cap Technology 
[17] for tuning Piccolo autopilot to be useful during the process of tuning the gains on 
Kestrel autopilot. The final PID gains set on Kestrel autopilot are presented in Appendix 
M. 


C; TEST SCENARIOS AND PROCEDURES 


Given the short duration of this thesis and this being the first attempt to run HIL 
for Kestrel Autopilot the experiments were conducted with the following decisions in 


mind. 


GPS waypoint navigation was used as the outer loop of the autopilot providing 


the guidance reference commands. 


The PID gains were obtained using all available knowledge the author has on 


Virtual Cockpit and Aviones. 

The aerodynamic performance parameters of the prototype were obtained from 
LinAir. 

Therefore, the results might not be a true representation of the expected 


performance of the prototype KUAV but the processes will be beneficial for future 


works. 
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1. Test Scenarios and Procedures 


A set of four waypoints (numbered 1 to 4) were set as the test route for the 
experiment. Waypoints | and 2 are set as intermediate positions for the UAV to reach the 
desired test altitude and speed. Waypoints 3 and 4 are on the actual test route, and 
waypoint 3 simulates the position where the KUAV is released from the HUAV. 
Waypoint 4 simulates the target position. Different release altitudes and horizontal ranges 
are achieved by adjusting the height and northing difference between waypoints 3 and 4 


respectively. 


A series of 7 runs, with velocities ranging from 6 m/s to 16 m/s at an interval of 2 
m/s, were conducted. (This range of test velocities is selected based on the results 
obtained from both the analytical and 6DOF + AP models in the previous chapters.) The 


detailed procedures are presented in Appendix N. 
D. RESULTS AND COMPARISON 


Using the procedures described in Appendix N, six simulations were completed 


and the result of a sample case (velocity of 8 m/s) is presented below. 
if Results 


Using the data generated from the HIL for a glide velocity of 8 m/s, the following 
graphs are plotted and shown in Figure 27. 


The graph of H vs. R, provides insights to the L/D, glide angle and maximum 
horizontal range traveled by the KUAV during its descend from 300 m. (Typical 


operating ceiling for tactical systems.) 


The graph of P vs. E, provides insights to the airborne endurance of the KUAV. It 
also shows the amount of transient the KUAU experiences after it is released from the 
HUAV. The average pitch angle will also be used to determine the AOA of the KUAV 
during its descend. The jump in pitch angle at the terminal phase is to be ignored as the 


KUAV had already landed on ground. 
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The graph of V vs. E, provides insights to the average velocity of the KUAV 
during its descend. It can also be observed from this graph that there is a transient period 
before steady-state condition is reached. The jump in velocity at the terminal phase is 


also to be ignored as the KUAV had already landed on ground. 
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Figure 27. HIL results for a gliding velocity of 8m/s 


The results for all the six simulations are summarized and attached in Appendix 


2 Comparison 


The results obtained from the HIL experiments were compared to those obtained 
from the analytical and the 6DOF + AP models. This comparison was only made for a 


small AOA range (until 4.5°) as this is the hardware’s physical limits during actual 


implementation. 
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Comparison Between Analytical, (DOF + AP and HIL Models for L/D vs. AOA 
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Figure 28. Comparison between analytical, (DOF + AP and HIL models for L/D vs. 
AOA 


Analysis of Figure 28 shows that the L/D computed from the HIL experiments 
displayed a similar trend but at a larger value as compared to those obtained in the 
analytical and 6DOF + AP models. In addition, the L/D value did not display a distinct 
point of maximum within this range of AOA. With the results from all 3 models 
displaying a similar trend, this further concludes the expected change of L/D with respect 


to changing AOA. 
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Comparison Between Analytical, (DOF + AP and HIL Models for R vs. AOA 
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Figure 29. Comparison between analytical, (DOF + AP and HIL models for R vs. 
AOA 


Analysis of Figure 29 shows that the R computed from the HIL experiments 
displayed a similar trend but at a larger value as compared to those obtained in the 
analytical and 6DOF + AP models. This result is expected, as it is shown in Figure 27 
that the values of L/D for HIL are larger. Similarly, the R value did not display a distinct 
point of maximum within this range of AOA. With the results from all 3 models 
displaying a similar trend, this further concludes the expected change of R with respect to 


changing AOA. 
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Comparison Between Analytical, (DOF + AP and HIL Models for E vs. AOA 
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Figure 30. Comparison between analytical, (DOF + AP and HIL models for E vs. 


AOA 


Analysis of Figure 30 shows that the E computed from the HIL experiments 
displayed a similar trend as compared to those obtained in the analytical and 6DOF + AP 
models. However, at higher AOA where the velocity of the glider is lower, the E obtained 
from HIL is significantly larger. This is clearly different from the expected decrease in E 


with increasing AOA as predicted from the analytical and 6DOF + AP models. 


The above results showed that the HIL model generates the same trend and 
behavior for EF and R with changing AOA as compared to the analytical and the 6DOF + 
AP models. However, the results obtained from the HIL experiments were better than 
expected. Ideally, the results obtained from the HIL experiments should fall in between 
the upper and lower bound predicted by the analytical and the 6DOF + AP models 
respectively. This shows that the simulation model used by Kestrel is more optimistic 
than of the analytical model. There is a need for future studies to better understand the 
control and guidance architecture used in Kestrel autopilot. Apart from the uncertainty 


mal 


caused by the autopilot, the aerodynamics parameters derived from LinAir might not be 
the actual representation of the KUAV’s performance. The performance data will have to 
be refined after physical flight test. A detailed list of recommended future studies is 


presented in the next chapter. 
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Vil. CONCLUSIONS 


A. CONCLUSIONS 


This thesis presented the concept of integrating COTS technologies to produce a 
prototype of a Tactical UCAV system. There is an existing HUAV developed and 
operated by the Naval Postgraduate School (NPS). For this reason, this thesis focused its 
effort on the design and evaluation of the performance of a non-powered glider with a 


wing that is linearly tapered and swept as the prototype KUAV. 


As there was no readily available model in the market that could fulfill the 
requirements of the proposed KUAV, a prototype was physically built for 
experimentation. The prototype was made purely from COTS hobby-grade material and 


components. The prototype was trimmed for low-altitude, hand-launched RC flight. 


The aerodynamic performances of the proposed prototype were first obtained 
using commercial panel code software. The aerodynamic performances were 
subsequently used in a Simulink 6DOF + AP software simulation, and the results 


obtained were comparable to those obtained from theoretical analysis. 


HIL simulations were conducted using the default setting in Kestrel Autopilot. 
This was the first attempt to run HIL testing for the Kestrel Autopilot. One objective of 
running the HIL was to validate the performance of the prototype prior to a high-altitude 
autonomous AP test flight. The other main objective was to establish the correct system 


architecture and procedures for future studies. 


Despite the short duration of this thesis, it shows that the proposed “Tactical 
UCAV” can be realized using COTS technologies. The following section describes some 


of the future studies equired to reach the envisaged end product. 
B. RECOMMENDATIONS AND FUTURE STUDIES 


The following future studies are recommended for developing the proposed 


“Tactical UCAV” system: 
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1. Understanding the Kestrel Autopilot 


There is a need to further understand the performance of the Kestrel autopilot. 
Given the short duration of this thesis, the current knowledge of the autopilot is 
superficial and insufficient to explain certain results or behaviors obtained from the HIL 


simulations. 
2. Variations to Gliding Profile 


This thesis had focused the research on investigating the maximum horizontal 
distance and airborne endurance the prototype can achieved given a released altitude. 
Future works should investigate the performance of the prototype for different glide 
profiles, for example the target is located right below the KUAV and not at a distance 


away. 
3. Physical Integration with ‘“‘Rascal’’ 


“Rascal” is an SUAV developed and operated by the Naval Postgraduate School 
(NPS), and it will be the HUAV carrying the KUAV. There is a need to develop a 
mechanical device to carry and release the prototype from “Rascal”. Only with this 


device can the prototype be tested with high-altitude launches. 
4. Data Integration with ‘Rascal’ 


Apart from developing the mechanical device to integrate the prototype with 
“Rascal,” there is also a need to establish an integrated method of communications for 
data transfers and updates between the HUAV and the KUAV. Two possible options are 


presented below for future works. 
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Hunter 


Option 1: Data transfers and updates are sent from the HUAV to the KUAV via 
Piccolo and Procerus ground stations as shown in Figure 31. 
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Figure 31. 


ground stations will be developed 


Communication protocol via ground control stations 


Option 2: Data transfers and updates are sent directly from the HUAV to the 
KUAV via XBee (Radio Frequency Module developed by MaxStream, Inc.) as shown in 


Figure 32. This option is preferred as the prototype KUAV need not have connectivity 









with the Procerus ground station. Effectively, the Procerus ground station can be 
Hunter 


eliminated and target data are relayed to the prototype KUAV through the Hunter. 


Target updates 
broadcast to PerceptiVU 









Continuous tracking of the 
target by the Hunter after 
the release of the Killer. 





Communication Protocol to 
link the hunter and the killer 
via XBee will be developed 













PVNT 


Workstation 











Killer 


TS 
’ 
Existing Comms 
Structure 



























The Killer glides to the 
\ 















































Target 
target based on the 
: updated target coordinates 
ATT Comput Vicia ——— 
empuer Piccolo Procerus ; 
ana NTS [**| Ground Station | | Ground Station | LOPtenal 
Ground Station ' 
Figure 32. 











Communication protocol via radio frequency module 


a) 


5s Variable Payload Capability for the KUAV 


The current aerodynamic performances of the prototype were obtained based on a 
fixed payload. The effects of variable payloads should be investigated, as the proposed 
KUAV is capable of carry variable payloads. 
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APPENDIX A. TECHNICAL SPECIFICATIONS OF PREDATOR 
B 





fal ed ry: h ge) 2B 


Persistent Multi-mission ISR 
and Strike Aircraft 








OBJECTIVE 
Perform multi-mission low- to high-altitude surveillance, 
reconnaissance, and “Hunter-Killer” missions. 


CHARACTERISTICS FEATURES 
Wing Span: 66 ft (20m) * Triple avionics 
Length: 36 ft (11m) * Redundant flight controls and control surfaces 
Powerplant: Honeywell TPE 331-10T * Remotely piloted or fully autonomous 
Max Gross Takeoff Weight: 10,500 Ib (4763 kg) * Military Standard 1760 Stores Management System 
Fuel Capacity: 4,000 Ib (1814 kg) * 6 wing stations for external carriage of payloads 
Payload Capacity: 850 Ib int. (385 kg) ¢ C-Band Line-of-Sight data link control 
3,000 Ib ext. (1361 kg) * Ku-Band Beyond Line-of-Sight/SATCOM datalink control 
Weapons: Hellfire missiles * Over 90% system operational availability 


GBU-12 laser-guided bombs —* Less than 12 hours displacement/emplacememnt 
GBU-38 Joint Direct Attack — * C-130 transportable (or self deploys) 


Munitions 
Power: 45 kVA, redundant 
Payloads: MTS-B EO/IR camera system 


Lynx SAR/GMTI 
Multi-mode maritime radar 


SIGINT/ESM system 
Communications relay 


PERFORMANCE 

Max Altitude: 50,000+ ft 
Max Endurance: 30+ hr 

Max Air Speed: 240+ KTAS 


me 
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APPENDIX B. 


TECHNICAL SPECIFICATIONS OF SPERWER B 





SPERWER B 
SPERWER: the European reference in TUAV 


Safe and airworthy design: 
+ Based on fail-safe and certified 
design 

+ Fully automated take-off and 
landing 

+ Advanced digital avionics suite 

+ Digital Ku band (15 GHz) data 
link 

+ Transponder/IFF mode 3C and 
VHF relay with ATC 


Tactical capacity: 

+ High Mobility 

+ Operations from unprepared 
areas 

+ Day and Night operations 

+ C4l connection and 
Interoperability 


Multi-payload capacity: 
+ EO/IR 

«SAR 

+ Transmission relay 

+ EW (ELINT, COMINT) 


Multi-mission capacity: 

+ ISTAR, targeting (20 m CEP) 

+ Long endurance (12+ hours) 

+ Weapon carrying capacity (60+ kg) 
+ Embedded sensor to shooter 
capability 


100% SPERWER 
compatible: 

+ Ground Segment compatibility 
- same Catapult pneumatic 
launcher 

- same Ground Control Station 
and Ground Data Terminal 

- same Maintenance concept 

- same Training 


+ Other UAVs compatibility 

- dual UAV capability from single 
system 

- combined operation with: 
SPERWER STD, Fast Penetrator 


59 





Overall length: 
35m 


Wing span: 

638m 

Ceiling: 

20,000 ft 

Speed: 

80 kts 

Data link range: 
200 km 

Payload capacity: 
100 kg 


re> Photos SAGEM, ETBS - MPC/D719 - 06/05 
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APPENDIX C. TECHNICAL SPECIFICATIONS OF SHADOW 200 





SHADOW 200: THE EYES OF THE COMMANDER 


AAI is a world leader in the design, production, and sustainment of tactical 
unmanned aircraft systems in support of warfighters and peacekeeping missions. 

Our Shadow 200s, known as the “eyes of the commander,” have exceeded 200,000 
flight hours, more than 85 percent in support of Operation Iraqi Freedom. Shadow's 
EO/IR payload sees targets 125 kilometers away from the brigade tactical operations 
center, providing situational awareness from an altitude up to 15,000 feet at more than 
3.5 kilometers slant range, day or night. 

Furthermore, our One System™ ground control station transmits imagery and teleme- 
try data directly to Joint STARS, the All Sources Analysis Systems, and the Army Field 
Artillery Targeting and Direction System In near real time, delivering an unequalled 
Interoperability network 
of intelligence-gathering 
capabllities. 

As a quick deployment 
asset, Shadow 200 tactical 
unmanned aircraft systems 
can be flown anywhere in 
theater using C-130 aircraft. 


































| U.S. Marines, who began transitioning 
to Shadow 200 systems in 2007, ready 
an aircraft for an intelligence-gathering 
sortie linked to a One System ground 
contro! station in the background. 





Ground crew can prepare a launcher and 
ready a Shadow 200 aircraft for flight in a 
matter of minutes. The hydraulic launcher 
mounts on a standard high mobility trailer. 





SYSTEM SPECIFICATIONS 


Performance Features Technical Data 

®@ Endurance: 5-7 hours ® Advanced Composite Construction @ Length: 11.2 ft. (3.41 m) 

@ Max Altitude: 15,000 ft. (4,572 m) ® Combat-Proven Hardware & Wingspan: 14 ft. (3.87 m) 

™ Max Dash Speed: 118 kts (218 kph) ® Rotary Engine with Proven Performance ™ Max Payload: 45 - 60 Ibs. (27 kg) 

® Cruise Speed: 90 kts (166 kph) © Multiple Payload Capability @ Max Gross Weight: 375 Ibs. (170 kg) 
® Loiter Speed: 60 kts (110 kph) © Multiple Launch and Recovery Methods 


§ Takeoff Distance (Launcher): 30 ft. (10m) ™ Operates with MOGAS 


® Modular Systern Design for Deployment 
Flexibility 
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APPENDIX D. TECHNICAL SPECIFICATIONS OF 
WATCHKEEPER TACTICAL UAV 


as 
Zins 


N 


SPECIFICATIONS - WATCHKEEPER TACTICAL UAV SYSTEM, UNITED KINGDOM 


DIMENSIONS 

Wingspan (WK180) 6.00rm 
Wingspan (WK450) 10.51m 
Length (WK180) 4.43m 
Length (WK450) 6.10m 
Fuselage Length (WVK180) 3.47mi 
WEIGHTS 

Take-Off Weight (WK180) 195kq 
Take-Off Weight (WK450) 450kg 
Payload (WK180) 35kq 
Payload (WWK450) 150kg 
PERFORMANCE 

Endurance (WK180) More than 10 hours 
Endurance (WK450) 20 or 30 hours 
Altitude (VK180) 4570m 
Altitude (WK450) 5 480m 
Loiter Speed (WK180) 130kméh 
Loiter Speed (WK450) 100km/h 
Maximum Speed (WK180) 194km/h 
Maximum Speed (WK450) 175km/h 
PROPULSION 

UEL AR 741 Engine (WK180) 28.3kVV 
UEL AR 801 Engine (WK450) 38.8kVV 
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APPENDIX E. COMPONENT SPECIFICATIONS 


A. RC Mode 





[S/No|Items | ‘Weight | Remarks 
a a gras | | 
RC Mode 
| i |2xSevo | | CCC 
| 2 |RadioReceiver | 6 | 
Battery Pack | 488 [ 
| 4 |PowerConverter | 6 | 
Nose Weight | 69 | lead 
Po Total} 77-7 [| 

Structural Weight 67.8 grams 

Total Weight 145.5 grams 

Effective Wing Area 71,438 mm? 


Average Wing Loading 2.04E-03 grams/mm? 
B. AP Mode 


[S/No|Items | «Weight | Remarks 
pT grams | 
| 1 |2xSevo | | 
| 2 lAutoPilot | 8 
| 3 |Modem | 88 | 


| 4 |GPSReceiver | 18.7 | 
| 5 |BatteryPack | 83-2 | 
| 6 |PowerConverter | 6 | 
[ee sGables =, «= ————— == ee 
| 8 |NoseWeight | 6.9 | Lead 
Pp Ctl? 106.6 | 





Structural Weight 67.8 grams 
Total Weight 174.4 grams 


Effective Wing Area 71,438 mm* 


Average Wing Loading 2.44E-03 grams/mm* 
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APPENDIX F. VERIFICATION OF PANEL CODES 


A. AR = | at Re = 70,000 
1. LinAir Input File (CD0 = 0.015 and K = 0.4) 


BD. AR1Glider - Notepad 
File Edit Format Yiew Help 
! Input data file for Linair 





! Reference values: 
Zref Mach Alpha #elements 
0.0 0.0162 1.90000 1 


! Element properties: 


! Element # 1 
Semi-Area  Semi-Span Sweep Dihedral 
10.0 31.2 0.0 
Root Inc. Tip Inc. 
0.0 .0 
# of panels 
10 





A Lindir Pro /-|o] x} 


File Edit Display Compute Help 
Reference Yalues 











Sref : 400.0000 
bref : 20.0000 
Href: 0.0000 
Yret: 0.0000 
Zref : 0.0000 
Nelen: 1 
alpha: 1.00000 
beta: g.00000 
phat: 0.00000 
qhat: o.00000 
rhat: 0.00000 
Mach: 0.01620 
WakeLoc: 1.000 
reflect: 1 
CLFile: NONE 
ForceFile: NONE 
ElementFile: NONE 
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Z Results generated by LinAir 


PB. AR1 - Notepad E ei) 


Edit Format View Help 





o 
i=] 


-00000 0.00000 
- 00000 00174 
00000 00347 
- 00000 00520 
00000 - 00694 
- 00000 - 00867 
00000 01041 
00000 -01214 
. 00000 - 01387 
00000 - 01560 
- 00000 - 01733 
00000 -01905 
- 00000 02078 
00000 -02250 
00000 + 02423 
00000 -02595 
00000 - 02766 
- 00000 -02938 
00000 - 03109 
- 00000 - 03280 
00000 - 03451 


lelejelejejeleleleselaelealelelelajeloajelas=) 
oo 
oo 
oo 
oo 
oo 


sielseielelelelelelelelelololeloses ola) a) 


oo 
oo 
oo 
oo 
oo 


SOOOOCOCOOCOOCCOCOOOCCOCOOD 
oo 
oo 
oo 
oo 
oo 


QO. 
QO. 
1. 
1. 
2. 
2. 
3. 
3. 
4. 
4. 
5. 
5. 
6. 
6. 
Pa 
va 
8. 
8. 
9. 
9. 
QO. 


CMOMOMOMoMouUdoMdUdUoOMdoD 
oooocoocoococococococoococ°c°c°c°c0oo 
oooococoococococococoococ°c°0c0o0ono 

Ooooococo0cooocooocooooooCoCoocCoo 
Ooooo0c0c0c0oocoocoooocoooocoo°coo 
lejesejejojojojejeaseslosesjojejesasjesoseal ala) 


1 














DOGS T T T T T T T T T 
I I I I i i \ 
\ I \ I I \ \ 
\ \ \ \ \ \ \ 
0.06}. — — — — — a jfSes= 1atoe= is SSeS \ | 
I \ \ \ \ \ \ \ 
\ \ \ \ \ \ I 1 \ 
\ \ \ \ \ \ 1 I \ 
OSS ae (enn el ae Atta, oes Ae aie, hie Py oes eae, | peer [Scien eked Waveney e en es, 
| I | I | I \ I Vv 
I I \ \ I | I \ J\ 
\ \ \ \ \ \ I 1 \ 
ooh — ~~ re Lp ee es cles ee eel eee ee eee (fees eee venienyeeeyeres| 
\ \ \ \ \ \ } ay, \ 
i 
/ 
ict onan ee ee Pi ie ha is oh ane ah ais le a ag Me a ae ae a Pe Mg teach 
fal \ 1 1 1 1 I \ J | 
s I \ \ \ \ \ \ a 
5 | | : : | | a | 
O08 Se Ieee oat a a OE fe es yo Ike ee ee ere 
‘9 \ \ \ \ \ \ YA 1 \ 
\ \ \ \ \ I 4 | i \ 
2 \ \ \ \ \ \ a \ 
[ed Fe ee [ee pee SS Pee Se aicasoa cea ae ieee [eS Sey fe BS Sy 
oO \ \ \ \ \ Wy 1 \ 
\ \ \ \ \ A I \ 
\ \ \ \ \ \ I \ \ 
OB eee [eee eos OR [ease ChE Se Deere [rene aires [pcb es eheeery execs iat nena eRe 
\ \ \ \ \ \ \ 1 \ 
\ \ \ | V \ 1 \ I 
\ 1 \ \ “| \ I \ \ 
0.025}- — — — — — In ----H joa -H 4qosf aH toon bos f nH ho---- In ----4 jos aH 4-----4 
\ \ \ \ 1 \ \ \ \ 
\ \ \ A \ \ \ \ 
\ \ \ \ \ \ \ \ \ 
oozh. — —— — — eee ere [eS dee Hens Pee eee Joo poe qe 
\ I i \ \ \ I I \ 
\ \ \ I \ \ \ \ \ 
\ i \ \ \ \ \ \ \ 
0.015 ain L i L L | L | | 
1 2 3 4 5 6 7 8 9 10 











AOA - Angle of Attack (deg) 


The above results, generated by LinAir for the AOA range of 0° to 10°, closely 
matched the results obtained in the wind tunnel experiments conducted by [12]. 
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AR = 2 at Re = 70,000 


I LinAir Input File (CDO = 0.020 and K = 0.4) 


P). AR2Glider - Notepad 
File Edit Format ‘Yiew Help 





! Input data file for Linair 


! Reference values: 
xref Mach Alpha #elements 
0.0 0.0 hs 0.0162 1.0 HE 


! Element properties: 


! Element # 1 
semi-Area  semi-Span Taper sweep Dihedral 
400.0 20.0 1 31.2 
xroot Yroot Zroot Root Inc. Tip Inc. 
0.0 0.0 0.0 0.0 0] 
cba cD1 cb2 # of panels 
0.020 i] 0.4 10 





Fj Lindir Pro /-|o) 
File Edit Display Compute Help 
Reference Yalues 


Sref : 800.0000 
bret : 40.0000 
Eref : 0.0000 
Yref : 0.0000 
Zret : 0.0000 
Nelemn: 1 
alpha: 1.00000 
beta: 0.00000 
phat: 0.00000 
ghat: 0.00000 
rhat: 0.00000 
Mach: 0.01620 
WakeLoc: 1.000 
reflect: 1 
CLFile: NONE 
ForceFile: NONE 
ElementFile: NONE 





69 


es Results generated by LinAir 


PB. AR2 - Notepad 
File Edit Format view 


Help 





[=] 
i=] 


SS0000000000C0COCDOCDDD 


QO. 
QO. 
1. 
HF 
2. 
2. 
3. 
3. 
4, 
4. 
5. 
5. 
6. 
6. 
7. 
7. 
8. 
8 

9. 
9. 
QO. 


DUHSMOUoUMoMnoMoMmounonouny 
SOooooooooCocCoooOoOoOoOoOooD 


1 


oooocooococoocoococoocoooceoo 


oooococococooocococcocooecocoeeoo 


0.00000 
0.02139 
0.04278 
0.06415 
0.08548 
0.10678 
0.12804 
0.14923 
0.17036 
0.19141 
0.21238 
0.23325 
0.25402 
0.27468 
0.29522 
0.31562 
0.33589 
0.35601 
0.37597 
0.39577 
0.41539 


0.02000 
0.02026 
0.02104 
0.02234 
0.02416 
0.02650 
0.02935 
0.03272 
0.03661 
0.04101 
0.04592 
0.05133 
0.05726 
0.06369 
0.07062 
0.07804 
0.08597 
0.09438 
0.10328 
0.11267 
0.12253 


0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 


0.00000 


-d. 
-a. 
-a. 
- 02903 
- 03484 
04065 
04647 
05229 
05811 
. 06393 
06976 
.07560 
. 08144 
. 08729 
- 09314 
09899 
-10486 
.11073 
.11661 


00580 
01161 
01741 
02322 


Q. 00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
o. 00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
Qo. 00000 
0.00000 


0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 


0.00000 
0.00361 
0.01389 
0.02941 
0.04830 
0.06872 
0.08918 
0.10868 
0.12659 
0.14268 
0.15687 
0.16924 
0.17996 
0.18918 
0.19708 
0.20382 
0.20957 
0.21444 
0.21855 
0.22199 
0.22486 





CD - Coefficient of Drag 
8 


2 
8 














5 


> 


OA - Angle of Attack (deg) 


The above results, generated by LinAir for the AOA range of 0° to 10°, closely 
matched the results obtained in the wind tunnel experiments conducted by [12]. 
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APPENDIX G. INPUT FILE FOR PANEL CODE 


f ginch_Glider - Notepad 
File Edit Format View Help 





! Input data file for Lindir 

! Reference values: 

larea Span xref yref Zref Mach Alpha #elements 
90516.0 501.8 0.0 0.0 0.0 0.0179 1.0 2 


Element properties: 


Element # 1 ; . 
Semi-Area  semi-Span Sweep Dihedral 


35719.0 250.9 7 34.3 0.0 
yroot Root Inc. Tip Inc. 
0.0 : 0.0 0.0 
cpl # of panels 
0.0 10 


Element # 2 
Semi-Area  sSemi-span sweep Dihedral 
il 29.9 0.0 
Root Inc. Tip Inc. 
0.0 0.0 
# of panels 
10 





Fj Lindir Pro | -(O} 


File Edit Display Compute Help 
Reference Yalues 


Sref: 90516.0000 
bref : 501.8000 
Kref : 0.0000 
Yref: 0.0000 
Zref: 0.0000 
Nelen: 2 
alpha: 1.00000 
beta: 0.00000 
phat: 0.00000 
qhat : 0.00000 
rhat: 0.00000 
Mach: 0.01790 
WakeLoc: 1.000 
reflect: 1 
CLFile: NONE 
ForceFile: NONE 
ElementFile: NONE 
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iy: 


APPENDIX H. DATA FILE GENERATED BY PANEL CODE 


Results for non-cambered airfoil 


BD. 8inf1 - Notepad 
Edit Format Yiew Help 


a 
a 





OVOWDOAINAMNMSPEWWNNEPR OO 


COSCO SCOSCOOCOOCOCOOCOOCOOCOOCOOCOOCOOCOOCOCOOCOCOOCOoCOoOCOSD 


DDODSCOOCOCOOOOOOOOOOCOOCOCOOCOOOCOOOOOOCOCOCOOOOCOODOOD 
DDODSCOOCOOOOOOOOOOOOCOCOCOOOOOOOOOOOCOCOCOOCOOCOOOODOOD 
DDODOSCOOCOOCOOOOOOOOOOCOCOCOOOOOOOODOOOCOCOCOOCOOCOODONDOD 
DODSCOOCOOCOOOOOOOOOOCOCOCOOOOOOOODOOCOCOCOOCOOCOOODODD 
DODD OCOCOCOCOOOOOOOOOOCOCOCOOOOOOOOOOCOCOOCOOCOOCOODONDD 
CODD OSOSCSCSCOOOCODCOOODOCOCOCOCOOCOCOOOODOOODOOOOCOOCOOOOODOOOOOD 
CODOSCOOCOCOOOOOOODOOOOCOCOOCOOOOODOOOOCOCOOCOOODOOOOODOD 


i} 
5 
0 
5 
0 
5 
a} 
5 
0 
5 
0 
5 
0 
5 
0 
5 
0 
5 
i} 
5 
0 
25 
-0 
25 
-0 
A] 
-0 
ae] 
-0 
25 
-0 
25 
-0 
25 
-0 
oe] 
-0 
a] 
.0 
23 
-0 
«3 
.0 


DOOSCOSCOSCOSCOCOSCOSCOSCOSCOCOCOCOCOOCOOCOCOOCOoCOOCOOD 
DOODOSCODCOSCOSCOOCOOCOOCOOCOCOOCOOCOOCOCOCOCOOCOCOOCOoCOoCOD 
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Results for cambered airfoil 


DB. Binf - Notepad 
File Edit Format View Help 


oO 
i=] 


0.05147 0.00000 -0.02463 
0.07718 0.00000 -0.03695 
0.10285 0.00000 -0.04928 
0.12848 0.00000 -0.06161 
0.15406 0.00000 -0.07396 
0.17958 0.00000 -0. 08632 
0.20501 0.00000 -0.09870 
0.23037 0.00000 -0.11109 
0.25562 0.00000 -0.12350 
0.28077 0.00000 -0.13594 
0.30580 0.00000 -0.14840 
0.33071 0.00000 -0.16089 
0.35547 0.00000 -0.17340 
0.38009 0.00000 -0.18595 
0.40455 0.00000 -0.19854 
0.42884 0.00000 -O.21115 
0.45295 0.00000 -O. 22381 
0.47688 0.00000 -O. 23651 
0.50061 0.00000 -O. 24925 
0.52413 0.00000 -0. 26203 
0.54743 0.00000 -0.27487 
0.57051 0.00000 -0.28775 
0.59335 0.00000 -0.30069 
0.61596 0.00000 -0. 31368 
0.63830 0.00000 -0.32673 
0.66039 0.00000 -0. 33983 
0.68221 0.00000 -0.35300 
0.70375 0.00000 -0. 36623 
0.72500 0.00000 -0.37953 
0.74596 0.00000 -0.39290 
0. 76662 0.00000 -0.40634 
0. 78696 0.00000 -0.41985 
0.80700 0.00000 -0.43344 
0.82670 0.00000 -0.44710 
0.84608 0.00000 -0.46085 
0.86512 0.00000 -0.47468 
0.88382 0.00000 -0.48859 
0.90217 0.00000 -0.50259 
0.92016 0.00000 -O. 51668 
0.93778 0.00000 -0. 53086 
0.95505 0.00000 -0.54513 


SCOODSCOSCOOSOOOSCOOSCOOCSCOSCOOCOODOOCSCOOCOCOOCOOOCOOCODD 


ecoocooocoocococoococococococoococococoococoocoooooooooooooooo 


OOWDANNAGMMbhWWNNRE ROO 


ie] 
5 
QO 
5 
9] 
5 
9] 
5 
8] 
5 
i¢] 
5 
is] 
5 
9] 
5 
9] 
5 
is] 
5 
8] 
a] 
-0 
EH 
.0 
23 
-0 
Ae 
0 
a] 
-0 
ae 
+0 
25 
-0 
a] 
-0 
23 
-0 
5 
-0 


OSOOOCOCOCOCOOCOOCOOCOCOCOCOCOOCOOCOOCOCOCOCOCOOCOCOOCOOoOoOoOoOoOoOoOoOo 
SCooococoococoococoocoooooooooooooooooooCoOOoOOoOoOOoOOoOoOoOoOoOo 
SCOODSCOSCOODOOCOODSCOOCSCOOCOOOCODOOOCOOCOOODOODODD 
eoocooocoococococococococococoocoocococococococococ0c0c0ococ°0c°oc°ocoo°0oo 
SOODSCOSCOOOOOCOOSCOSCSCOOCOOOCOOCOOSCOOOOOCOOOCODOD 
lejelelelelelejejejesejesejejejejelejesesesjesesesesejeseajejesejesejesesesejeleleal =) 





Column | = Angle of Attack 

Column 2 = Angle of Sideslip 

Column 3 = Free Stream Mach Number 
Column 4 = Lift Coefficient 

Column 5 = Drag Coefficient 

Column 6 = Sideforce Coefficient 
Column 7 = Moment Coefficient 
Column 8 = Rolling Moment Coefficient 
Column 9 = Yawing Moment Coefficient 
Column 10 = Span Efficiency 
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APPENDIX I. = INITIALIZATION FILE FOR 6DOF + AP MODEL 


fo) 


& Initialization file for 6DOF + AP simulation model. 





eCLe: 
T=0.1; 
r_lim=7; 


& Initial Conditions in ENU (all vector data is represented as a column 
vectors) 











Pos_0 = [0; 0; 300]'; % Initial position vector (m) 
Buler_O = [0; 0; O]'*pi/180; 6 Initial Euler angles (rad) 
Omega_O = [0; 0; O]'; & Initial Omega (rad/s) 
POQR_O = [0% 0: OQ] '; % Initial Omega (rad/s) 
Vb_0O = [5; 0; O]'; % Initial body-velocity vector (m/s) 


% Mass and Geometric Parameters recomputation 





S (0.2032+0.1576) *0.5018/2; % Surface area of wing (m2) 
span 0.5018; %& wingspan (m) 
chord = 0.1814; % Mean Aerodynamic Chord (m) 
mass = 0.175; % gross weight (kg) 
Ixx = 0.0028959; %& main moment of inertia around axis Ox (kg*sq.m) 
Iyy = 0.0006514; & main moment of inertia around axis Oy (kg*sq.m) 
1z2z = 0.0035430; %& main moment of inertia around axis Oz (kg*sq.m) 


% Aerodynamic Derivatives (all per radian) 

CLO = 0.05147; lift coefficient at a = 0 

CLa = 2.9095; lift curve slope 

CLa_dot = 0; lift due to angle of attack rate 
Cc 

Cc 

Cc 


ae ol? 


ole 


im 
QO 
ll 

© 
~ 

oe 


lift due to pitch rate 
lift due to elevator 
drag coefficient at a = 0 





C 
0 
oO 
| 
| 
at 
W 
~ 
~ 
Oy 
~e 
oe 





> 
oO Oo 
Oo 
ae 

ie) 

6 

od 

OO 
oe) 
N Bd 
Nh Nb 
OO 
weoN 
ae ol? 





























drag curve slope (A2) CD=CD0+A2*1.25*CL*2 
Al = 0; 
CYb = 0; % Side force due to sideslip 
CYDr = 0; & sideforce due to rudder 
Clb = -0.00172; % dihedral effect = -0.00172 
Clp = 0; % roll damping 
Clie = 0; % roll due to yaw rate 
Cl1Da = 0; % roll control power 
ClDr = 0; % roll due to rudder 
Cm0 = -0.02463; & pitch moment at a = 0 
Cma = -1.4192; %& pitch moment due to angle of attack 
Cma_dot = 0; % pitch moment due to angle of attack rate 
Cmq = 0; % pitch moment due to pitch rate 
CmDe = -1.1543; % pitch control power 
Cnb = 0; % weathercock stability 
Cnp = 0; % adverse yaw 
Cnr = 0; % yaw damping 
CnDa = 0; % aileron adverse yaw 
CnDr = 0; % yaw control power 
CLD£E = -1.3776; % lift due to flap 
CmDf£ = -1.1543; % flap control power 


7S 


fe) 


%* Standard Atmosphere 





ISA_lapse = .0065; % Lapse rate (degC /m) 
ISA_hmax = 2000; % Altitude limit (m) 

ISA_R = 287; % Gas Constant (degk*m*m/s/s) 
ISA_g = 9.815; %& Gravity (m/s/s) 
ISA_rho0O = 1.225; % Density at sea level (kg/m/m/m) 
ISA_PO = 101325; % Sea-level Pressure (N/m/m) 

ISA_TO = 289; % Sea-level Temperature (degK) 
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APPENDIX J. SIMULATION RESULTS FROM THE 6DOF + AP 
MODEL 


[8.0 | 105.0 
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APPENDIX K. PROCEDURES FOR SETTING UP THE TCP/IP 
CONNECTION 


The following set of procedures is used to set up the TCP/IP connection for the 


HIL experiments. 
Remove the modem from the autopilot. 


Disconnect the programming wire on the programming cable as shown in Figure 


a3 





Figure 33. Disconnected programming wire 


Plug the 5 pin header of the programming cable into the autopilot modem port as 
shown in Figure 34. The top cable is the programming cable, and the bottom cable is the 


power cable. 





Figure 34. Connection for the programming and power cables on autopilot 
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Plug the other end of the programming cable into a free serial port on the 


computer that will be running Virtual Cockpit. 
Power on the autopilot. 
Load Aviones. 
Load Virtual Cockpit. 


Open the “Edit Agents” window in Virtual Cockpit. (This can be found in the 
“Agent Menu”) The right-most column is labeled “HIL” for hardware-in-the-loop. 


Check the box next to the agent number of the autopilot as shown in Figure 35. 


Edit Agent List 


on | Address Name} Color HIL| st| 
lv 1032 AVI Iv | 
[ Hii 





Figure 35. HIL selection in Virtual Cockpit 


Click “OK.” 


A dialog box will pop up (see Figure 36) in which the IP address of the computer 
running Aviones may be specified. If Aviones is running on the same computer that is 
running Virtual Cockpit, click the “Loopback” button and click “OK” to connect to 
Aviones. Virtual Cockpit will now attempt to connect to Aviones and add all the selected 


agents that had the HIL checkboxes checked. 


Connect to Avionies 


IP Address 
1228 











Figure 36. Connection window for Aviones 
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If connected successfully, the Virtual Cockpit will indicate on its status bar and 
Aviones will have added the Agent to its list on the bottom of the window, as shown in 


Figure 37. 


@@ Aviones SEE 


File UAV View Physics Control Help 


r= 


= 





Agent Time Easting | Northing | Alt HAG Airspeed | Heading 
5.65 oO Q 50 50 10.0 0.00 


Figure 37. Agent added in Aviones 


The next step is to make sure that Aviones is set up to talk to the autopilot through 
the correct comm port. Open up the HIL simulation control dialog in Aviones through 


“View Menu” followed by “HIL Sim Ctrl’. 
Place the cursor in the Serial Port # box. 


Delete the value that is currently in the box using the backspace key. (The box 
should turn red) 


Enter the serial port number connected to the autopilot. 


Click “Enter” — The box should turn white and the word “Open” should be 
displayed in the box to the right as shown in Figure 38. 
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HIL Simulation Control [x] 


Agent Number | Agent1032 A 


Serial Information 


Serial Port # fi Open 
Open | Close | 


Starting Latitude Starting Longitude 
0,000000 0,000000 


Starting Altitude 


50,000000 


Commands 


Loss of GPS | Loss of Comm 
Engine Failure | 


Figure 38. HIL simulation control window 





With the correct serial port open, data should begin to flow between Aviones, 


Virtual Cockpit and the autopilot. Virtual Cockpit should show communications with the 


autopilot. 


To verify communication between Aviones and the autopilot, the right-most LED 
on the autopilot should be blinking rapidly. This indicates that the autopilot is receiving 


sensor information from Aviones. 


The final display on the user interfaces for Aviones and Virtual Cockpit (if there 


is a successful connection) is shown in Figure 39. 
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Ga [ant (es [HE] [we 


tne} ED (ns) (om 


‘Address| Comm Batt | Akinude | Velocty | Sat] Rc 





Figure 39. User interfaces for Aviones and Virtual Cockpit 
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APPENDIX L. 


PROCEDURES FOR UPDATING AERODYNAMIC 


PERFORMANCE PARAMETERS IN AVIONES 


The initial performance parameters loaded in Avoines were for Zagi UAV test 


platform. The procedures for updating these performance parameters are as follows. 


Locate the initialization text file, “physics params” in the Aviones folder. 


Change the necessary parameters to match the performance of the prototype. The 


updated parameters for this thesis are shown in Figure 40. 


Figure 40. 


P physics_params - Notepad 
File Edit Format View Help 


jG) 





gravity, 9.806650 
5 


$,0.12945 

S prop, 0.033400 
cxdt,1. 000000 
Kmotor, 25. 000000 

cL 0,0.051470 

cD 0,0.022000 

cm 0,-0.02463 

CL alpha, 2.909500 
CD alpha, 0.688700 
cm alpha, -1.4192 

cL gq, 0.000000 

cD q, 0.000000 

cm q,-3.600000 

CL elevator, -1.377600 
cD elevator,o. 000000 
cm elevator, -1.154300 
cy 0,0. 000000 

cl 0,0. 000000 

cn 0,0. 000000 

cY beta, -0. 980000 
Cl beta, -0.020000 
cn beta, oO. 250000 

cy p,o. 000000 

cl p,-0. 200000 

cn p, 0.022000 

cy r,0. 000000 

cl r,0.140000 

cn r,-0.350000 

cy aileron, o. 000000 
cl aileron, o. 080000 
cn aileron, o. 060000 
cy rudder, -0.170000 
cl rudder, 0.105000 
cn rudder, -0.032000 





(>| 


li 





Aerodynamic performance coefficients for prototype 
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Save the changes. (Do not change the file name) 
Run Aviones. 


Verify that the parameters had been updated by clicking on “View” followed by 
“Physics Parameters” as shown in Figure 41. The values of the parameters can also be 


changed from this user interface. 


M@ Aviones 


| File UAV | Niet Physics Control Help 
UAV 1 
UAV 2 
UAV 3 
UAV 4 
UAYS 
UAV 6 


Physics Parameter Editor x 





Value 
9.806650 
0.175000 
0.002896 
0.000651 
0.003543 
0.001500 
1.422400 
0.181400 
1.268200 
0.129450 
0.033400 
1,000000 
25 AnnANn 








Chase Perspective 
Onboard Camera Perspective 


Satellite Perspective 
View UAY from Launch Pos 
View All UAYs From UAY 1 





Plot Window 
Plot Window Options... 








Autopilot Registers... 
Physics Para 
Control Window... 





HAG Airspeed | Heading 





HIL Sim Control 





Figure 41. Physics parameter editor window 
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APPENDIX M. PID GAINS FOR PROTOTYPE ON KESTREL 


AUTOPILOT 


The PID gains can be changed from the user interface of Virtual Cockpit. Click 


on “Settings” followed by “PID Values”. 


The PID gains for the prototype on Kestrel autopilot are shown in Figure 42. Only 


the parameters shaded in yellow was changed. The remaining parameters were the default 


values used for Zagi UAV test platforms. 


G Global Settings 





UAV-1032 (F5) 


Autopilot Config 
Autopilot Variables 
Data Logs 
Failsates 

» PID Values 
Autopilot Sensors 
Servos 
Calibration 
Payload 

Gimbal 

System 


UAV-1032 (F5) 


a 





Level 3 Loops Level 2 Loops Level 1 Loops 








Roll Desired 





Actual | 





Hdg Fixed Input Foll Fixed Input 
Pitch Dyn Input 


O Altitude 
© Airspeed 
[J Auto Alt 





0,000 ° 
15.986 ° 
0,000 ° 
6,000 m/s 
100.000 m 
0,000 rad/s 
0,000 m/s 


Roll 
Pitch 
Heading 


Roll Rate 
Pitch 
Pitch Rate 
Yaw Rate 


O Thr>Airspd 
O Thr>Alt 
© Thr->Climb/Alt 














Kd 
Vv 
v 
(VJ 


Tuning 

(| Use Desired - 
(Disable 
Joy/Kepbrd 
Packets] 

















Airspeed 
Altitude 
Turn Rate 
Climb Rate 











0,000 ° 
0,000 ° 
0.000 ° 
9,950 mis 
50.000 m 
0,000 rad/s 
-0,063 m/s 





Y-Max |-1 Y-Min [¥] Auto-X [ 








Upload Loops||Update Flash) { foe | (Cas | 1 


Auto-Y 











Pitch Kp 

Pitch Ki 

Pitch Limit 

Pitch Rate Kp 

Pitch Rate Limit 

Roll Kp 

Roll Kd 

Roll Ki 

Roll Limit 

Roll Rate Kp 

Roll Rate Limit 

Throttle <-Airspd Kp 
Throttle<-Airspd Kd = 1 
Throttle<-Airspd Ki 0.1 
Throttle <-Airspd Limit 100 % 
0.15 
Yaw Rate Limit 0.35 rad 


Yaw Rate Kp 





Thr<-Climb RateKp 3 


Figure 42. Values of PID gains 
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APPENDIX N. PROCEDURES FOR RUNNING THE HIL 
EXPERIMENTS 


As mentioned in Chapter VI, this was the first attempt to run HIL for Kestrel 
Autopilot. Therefore, this set of procedures is described in detail for the benefit of future 


works. 
Set up the TCP/IP connection as described in Appendix J. 


Create a test route for the respective velocity (6 m/s to 16 m/s at 2 m/s intervals), 
with waypoints 3 and 4 having a height difference of 300 m. Upload the test route to the 
autopilot using the flight plan control window (click on the red “Upld” tab) of Virtual 
Cockpit as shown in Figure 43. 


fa (=) [ios] (=) [a 


(ems) [ Nev 


Address | Coenen | 





Figure 43. Flight control window of Virtual Cockpit 
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Launch the UAV using the “UAV” tab in Aviones by clicking on “Launch if 


necessary” or “Relaunch current UAV” as shown in Figure 44. 


@ Aviones 
File | EW View Physics Control Help 
Add new... 





Relaunch current UAY 
Instantiate UAV and launch 


Time Easting | Northing | Alt HAG Airspeed | Heading 
264,39 o 0 50 50 10.0 0.00 








Figure 44. Launching the UAV from Aviones 


Ensure that the UAV is in “Nav” mode and the autopilot is navigating to 


waypoint | as shown in Figure 45. 
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aProcerws Virtual Cockplt (~ 51x] 





ia lacie 


fr [Maru] | 


Address | Comes Wart Atul Velocky | Sat| RC 
msl 6 





Toot Fite tpt MQ GG} (Deven) (Cold) Cine wl» 





Page rene hee be eee Ure rise [rs] Takeo Sod (m/s) Alt (mn) Spiral Radius North East Latitude Longitude 
0.00. He mn: wm cd 


— | Sam: Connected, Dev: No Connect KEYBOARD 


Figure 45. UAV in“NAV’ mode 


The benefit of using Virtual Cockpit’s default settings for the Zagi UAV test 
platform is that the Zagi is a powered glider. It is capable of autonomous flight toward 


the active waypoint with its onboard propeller motor. 


When it has reach waypoint 3, which simulates the release point of the KUAV, it 
is necessary to “turn off’ the propeller motor as the KUAV is a non-powered glider. This 
can be implemented through Aviones. Under the “View” menu in Aviones, open “HIL 
Sim Control” and click on the “Engine Failure” tab as shown in Figure 46. This will 


“tum off’ the motor, but the autopilot is still navigating waypoints. 
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@® Aviones Ce) HIL Simulation Control (X) 






Agent1032 (1032) Agent Number Jagent1 032 7 | 
UAV 2 + Serial Information} 
UAV 3 SerialPort # [1 [Open 

UAY 4 

UAVS Open | Close | 











UAV 6 

















Chase Perspective Starting Latitude Starting Longitude 






Onboard Camera Perspective Jo.oon000 Jo.oog000 
Satelite Perspective 0.000000 0.000000 
View UAV From Launch Pos 


Starting Altitude 


| 50,000000 


Loss of Comm. | 


Low Battery 






View All UAVs from UAY 1 


Plot Window 
Plot Window Options... 


Autopilot Registers... 
Physics Parameters... 
Control Window... 







Commands 
Engine Failure 





















HIL Sim Control 


Figure 46.  HIL Sim Control Window of Aviones 


Virtual Cockpit has a data logging capability and this is extremely useful in 
recording the results for subsequent analysis. The data log can be activated by clicking 


“Settings” followed by “Data Logs”. The user interface is shown in Figure 47. 


G Global Settings 






























Precision Data Logger 
























































UAV-1032 Enter search here... | Clear Logged Variable List ice eaeCE 
Settings Square Effort Scalar a Airspeed 7 =| Sa 
Autopilot Config Status to GCS Send Time az Altitude 2 100 Samples 
‘Aulopaot Variables System Flags? GPS Notth Position = 
Pe System Flags2 Theta 20 Seconds 
» Data Logs System Flags3 
Failsafes System Status] 
Take Off STage 
PID Vaiss Take Off Timer 
Autopilot Sensors Temp Comp Mode 
Servos Temperature R Download Log: 
Calibration Terrain Elevation 100 bytes at a time 
Theta Quat 
Payload Throttle Servo Position E AHI MO Ae 
Gimbal Time Since Last Comm 5 — — 
System Time Since Last Good Battery | 
Time Since Last Good RPM 
Time Since Last GPS Comm Download 
Time Since Last GPS Lock 
Time Since Last RC Over Comm 
Time Since Last RPM Comm 
Time Since Launch 
True Airspeed 
True Psi 
UAV Mode 
Update GPS 
‘Wind Heading 
‘Wind Speed 
Yaw Gyro rv Save. 
rape Dorcel anr 
Telemetry Logger 
Base Filename:  telemetry_%n_%Y_%m_%d_#.m M-file (MatLab Plain Text) 
telemetry_UAV1_2007_11_19_00.m 1376 bytes 0 samples - ck 
Start 
UAV-1032 (F5) = 











Debug Packet 


Cig 


Figure 47. Data log window of Virtual Cockpit 
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Data logging commences moments before engine failure is activated to simulation 


the releasing of the KUAV from the HUAV. 
Four data are recorded for analysis. 


Airspeed, the speed at which the KUAV is descending after being released from 
the HUAV. 


Altitude, the height of the KUAV during its descend after being released from the 
HUAV. 


GPS North Position, the horizontal range which the KUAV had traveled during its 
descend after being released from the HUAV. 


Theta, the pitch angle of the KUAV during its descend after being released from 
the HUAV. 


When then altitude of the UAV reaches zero, this indicate that the UAV have 
landed. 


The data for each run are saved as M-scripts and analyzed using Matlab. The 


following script is written to analyze the results and plot the respective graphs. 


sData analysis for HIL 
cle; 
test8ms; 


data=datalog.samples; 


E=data(:,1)./5; 
P=data(:,3).*180/pi; 
R=data(:,4); 
V=data(:,5); 
H=data(:,6); 





glide=(data(100, 6) -data(250,6))/(data(100, 4) -data (250, 4) ) 
glide_deg=glide*180/pi 





pitch=mean (data (100:250,3) ) 
pitch_deg=pitch*180/pi 


aoa_deg=pitch_deg-glide_deg 


Velocity=mean (data(100:250,5) ) 
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Start_Time=data(1,1); 
End_Time=data(316,1); 
Endurance=(End_Time-Start_Time) /5 





Start_R=data(1,4); 
End_R=data (316,4)j; 
Range=(End_R-Start_R) 





figure (1); 
hold on; 





subplot (3,1,1) 








plot (R,H); 

axis equal; 

xlabel('R - Horizontal Range (m)', 'fontsize',14); 
ylabel('H - Altitude (m)','fontsize',14); 








title ("H ws. R', 'fontsize',16¢); 
grid on; 


subplot (3,1,2) 

















plot (E,P); 
xlabel('E -—- Endurance (sec) ','fontsize',14); 
ylabel('P - Pitch Angle (deg)','fontsize',14); 

















title('P vs. E', 'fontsize',16); 
grid on; 


subplot (3,1,3) 























plot (E,V); 
xlabel('E -—- Endurance (sec)', 'fontsize',14); 
ylabel('V - Velocity (m/s)','fontsize',14); 











title("¥ we, EB", 'fontsige*, 16); 
grid on; 
print -tiff -depsc 8ms.eps; 


hold off; 


Repeat the above procedures for the other test velocities. 
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APPENDIX O. SIMULATION RESULTS FROM HIL 


Saas EE a Ee a En eS eS es 
|__ (deg) | (deg) | ~~ (deg) — | (mT (s) Svs), Tm) 
| 5.06 | 9.75 | OT t752 | 8 OT 00S 584 
| 8.47 10.99 | 252594 | 200 | 8800S 5.31 
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